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g : ABSTRACT 
O 

We present the first determination of the 15//m luminosity function of galaxies 
from the European Large Area ISO survey (ELAIS) southern fields. We have 
adopted a new criterion to separate the quiescent, non-evolving and the starburst, 
evolving populations based on the ratio of mid-infrared to optical luminosities. 
Strong evolution is suggested by our data for the starburst galaxy population, 
while normal spiral galaxies are consistent with no evolution. The starburst 
{Sj | population must evolve both in luminosity and in density with rates of the order 

^ ; L(z) oc (1 + z) 3 - 5 and p(z) oc (1 + z) 3 - 8 up to z ~ 1. The evolutionary parameters 

of our model have been tested by comparing the model predictions with other 
observables, like source counts at all flux density levels (from 0.1 to 300 mJy) 
and redshift distributions and luminosity functions at high- z (0.7 < z < 1.0 from 
i- £h HDF-N data). The agreement between our model predictions and the observed 

data is remarkably good. We use our data to estimate the star-formation density 
of the Universe up to z=0.4 and we use the luminosity function model to predict 



the trend of the star-formation history up to z — 1. 
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1. Introduction 

The extragalactic background light shows that the emission from galaxies at infrared 
and sub-millimeter wavelengths is an energetically significant component of the Universe. 
This emission originates from star-formation activity and active galactic nuclei. The precise 
contribution from each type of activity is still debated. It is thus important to our under- 
standing of galaxy and AGN formation to study those populations that emit a substantial 
amount of light at infrared (IR) wavelengths in their rest frame. In particular, data from 
deep Infrared Space Observatory (ISO) surveys at 15 /im (i.e. Elbaz et al. 1999; Flores et al. 
1999; Lari et al. 2001; Metcalfe et al. 2003) seem to require strong evolution for galaxies emit- 
ting in the infrared wavebands. This result, supported also by the detection of a substantial 
cosmic Infrared Background in the 140 /im - 1 mm range (Puget et al. 1996; Hauser et al. 
1998; Lagache et al. 1999), has stimulated the development of several evolutionary models 
for IR galaxies (i.e. Rowan-Robinson 2001, Franceschini et al. 2001, Chary & Elbaz 2001, 
Xu et al. 2003). All these models fit with different degrees of success the IR/sub-millimeter 
source counts and the cosmic Infrared Background within the present uncertainty limits, 
but suffer of parameter degeneracy and none of them is based on a local luminosity function 
(LLF) obtained from 15-/im data, being all extrapolated from different IR wavelengths (12, 
25, 60 /im). 

So far, complete spectroscopic samples of 15-/im sources have been obtained only in 
small fields (i.e. HDF-N: Aussel et al. 1999; HDF-S: Mann et al. 2002, Franceschini et al. 
2003), too small and too deep to allow a detailed study of the local luminosity function. 

The ELAIS survey is the largest open time project conducted by ISO (Oliver et al. 
2000), mapping an area of ph 12 deg 2 at 15 /im and 90 /im. The final, band- merged ELAIS 
catalogue has recently been completed (Rowan- Robinson et al. 2003). The spectroscopic 
data is most complete in the southern fields and in this paper we present an analysis of the 
15-/im luminosity function derived from these data. This is the first determination of the 15- 
/im luminosity function and its evolution, constrained by all the available observables in this 
band (source counts from IRAS to the deepest ISOCAM flux densities; redshift distributions 
at low and high- ,2 using both data from ELAIS and data from the deeper HDF-N survey). 
The model fitting the LF is then used to estimate the star-formation history of the Universe. 

This paper is structured as follows. In Section 2 we present our data sample. In Section 
3 we discuss the adopted IR and optical K-corrections. In Section 4 we show the method 
used to compute the 15-/im LF and present the results. In Section 5 we compare our LLF 
determination with previous ones. In Section 6 we discuss the evolution rates derived from 
our data, and compare the model predictions with the observable constraints and with other 
literature models. In Section 7 we present our conclusions. 
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Throughout this paper we will assume H = 75 km s 1 Mpc x , Vl m = 0.3 and = 0.7, 
unless explicitly stated. 

2. The data samples 

2.1. The parent catalogues 

Our analysis uses the southern ELAIS fields, which have more complete spectroscopy 
than the northern fields at the present time. SI and S2 are the main survey fields in the 
southern hemisphere: SI is centered at a(2000) = 00 A 34 m 44.4 s , 5(2000) = -43°28'12" 
and covers an area of 2°x2°, while S2 is centered at a(2000) = 05 ft 02 m 24.5 s , 6(2000) = 
-30°36'00" and covers an area of 21'x21'. 

The 15-/zm data in these fields have been reduced using the LARI technique described 
in detail in Lari et al. (2001) and the 15-/zm catalogues obtained in the two fields are given 
by Lari et al. (2001) and Pozzi et al. (2003), respectively. The source counts derived from the 
main field SI are presented and discussed in Gruppioni et al. (2002). The whole SI and S2 
areas have been surveyed in the radio down to Si AGHz w 0.2 and 0.13 mJy respectively. The 
radio data analysis is presented in Gruppioni et al. (1999) and Ciliegi et al. (in preparation). 

SI contains 462 sources detected at 15-/zm, 406 of which constitute a highly reliable 
sample which can be used in a conservative statistical analysis (see La Franca et al. 2004), 
while S2 contains 43 15-/iin sources. Optical photometric follow-up has been obtained in 
SI with the ESO/Danish 1.5m Telescope (to R ~ 23) and in S2 with the ESO WFI/2.2m 
Telescope in U, B and / (to / ~22) and SOFI/NTT in K' (to K' ~ 18.75). The spectroscopic 
follow-up was carried out at the 2dF/AAT and ESO 3.6m, NTT and Danish 1.5 Telescopes. 
In SI, ~81 % (328/406) of the 15-/im sources have an optical counterpart in the R band 
(taking into account all the optical identification, sources to R ~23 plus 6 sources with 23 
< R < 24.3, see La Franca et al. 2004), while in S2, where the identification has been done 
in the I band, the percentage is higher (~90 %, 39/43). In SI, the spectroscopic redshift 
completeness is ~88 % of identified objects and ~71 % of the whole sample (290 objects), 
while in S2 the corresponding figures are ~77 % and ~70 % (30 objects). 

2.2. The spectroscopic sample 

To study the 15-/zm luminosity function and its evolution we restrict our analysis to 15- 
/iin sources with optical counterparts in the ranges of magnitude with higher spectroscopic 
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identification completeness. Recalling the definition given in La Franca et al. (2004) of the 
Sljrest region (the shallower area) and of the Sl_5 region (the deeper part), we consider the 
following objects: in Sl_rest all the sources down to R < 20.5 (97 % complete, 216/223); 
in Sl_5 the sources down to R < 21.6 (97 % complete, 73/75); in S2 the sources down to 
R < 22.6 and with 15-fim fluxes >1 mJy (100 % complete, 18/18). In S2, the magnitude 
limit of / ~ 22 can be converted (for uniformity with SI) to a limit in R magnitude, by 
assuming a mean colour R — I ~ 0.6 (found for our galaxies). This leads to an R limit of 
~ 22.6 (see Pozzi et al. 2003). The conservative 15-/im cut is adopted in S2 because in this 
field there are uncertainties in the 15-/im completeness function below this flux level. 

The sample of galaxies with measured redshift, after excluding objects classified as AGNs 
(both type 1 and type 2) or stars on the basis of their spectra, consists of 161 galaxies: 
101 in Sl_rest, 48 in Sl_5 and 12 in S2. Following La Franca et al. (2004) and Pozzi 
et al. (2003), AGNs have been separated from galaxies according to classical diagnostic 
diagrams and information on the nature of galaxy objects has been provided following the 
classification scheme of Dressier et al. (1999). This scheme is based primarily on two lines, 
[OII]A3727 in emission and H<5 in absorption, which are good indicator of (respectively) 
current and recent star formation. Our sample of 161 galaxies include 40 e(a) (spectra of 
dust-enshrouded starburst galaxies), 13 e(b) (spectra with very strong emission lines), 72 
e(c) (spectra typical of spirals), 32 k(e) (spectra with signs of at least one emission line) 
and 4 k (spectra of elliptical-like objects) galaxies. In Figure 1 we show the magnitude 
vs. 15-fim flux distribution of the whole Sl_rest+Sl_5+S2 sample of galaxies (objects with 
spectroscopic identification are shown with filled symbols). 

In our estimates of the luminosity function and its evolution (Section 4) we apply an 
additional restriction, selecting only galaxies with redshift z <0.4. As shown in Figure 4, in 
fact, above this redshift threshold the redshift distribution of galaxies is poorly sampled. 
The final total sample of galaxies considered for our LF determination consists of 150 objects. 

3. The K-correction 

To derive the 15-/zm luminosity function from our sample of galaxies we need to relate 
the observed flux S (and redshift) to the rest-frame 15-/xm and i?-band luminosities. This 
requires the knowledge of the K-correction in the Mid-IR and in the i?-band. 

Models of infrared emitting populations (Rowan-Robinson 2001; Franceschini et al. 2001; 
Elbaz et al. 2002; Xu et al. 2003) have shown that the IR counts can be explained with three 
main populations. These populations are characterised by their SED and their evolutionary 
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properties: the normal spiral galaxies, the starburst galaxies and galaxies powered by an 
active galactic nucleus (AGNs). The three populations can be separated on the basis of 
their SED in the Mid-IR/Far-IR band (see Rowan-Robinson & Crawford 1989; Xu et al. 
2003). 

We have assumed M51 and M82 as prototypes of the normal and the starburst pop- 
ulations respectively (the work on AGNs is presented in Matute et al. 2002; Matute et al. 
in preparation). M82 is a local moderate-starburst galaxy (Ljr ~ 1O 1O ' 6 L , Elbaz et al. 
2002). We did not consider the SEDs of a more active galaxy (like Arp220, L//j^10 11 Lq), 
since from the identifications of ELAIS (La Franca et al. 2004) and HDF-N sources (Elbaz 
et al. 2002), such extreme galaxies are rare and tend to appear only at high redshift (z ^0.8). 
The model SEDs of both M51 and M82 have been taken from the GRASIL model output 
(Silva et al. 1998) except for the Mid-IR (5-18 fxm) spectrum of M82, for which the observed 
ISOCAM CVF (Circular Variable Filter) spectrum (Forster Schreiber et al. 2003) has been 
adopted. The two reference SEDs (normalized to the bolometric luminosity from 0.1 to 1000 
jum) and the corresponding K-corrections in the Mid-IR filter (ISOCAM LW3 filter) and in 
the i?-band filter (R Cousin) are shown in Figure 2. The normal spiral galaxy has relatively 
high Far-IR over Mid-IR emission, and has a bolometric luminosity dominated by the optical 
region of the spectrum. The active galaxy has a bolometric luminosity dominated by the 
IR emission, as expected in presence of star-formation obscured by larger amount of dust. 
As shown in Figure 2 (bottom), the two galaxies have similar K-corrections in the LW3 and 
i?-band filters, at least up to z~l. 

Since in our sample we do not have the Mid/Far-IR colours, we have used the Mid- 
IR to optical luminosity ratio to tentatively associate a representative SED (M51 or M82) 
to each galaxy. As discussed by La Franca et al. (2004), galaxies selected in the Mid-IR 
band show a well defined relation between the ratio of mid-infrared to optical luminosities 
(Li 5 ^ m /L R ) and the mid-infrared luminosities. This relation, although with some significant 
scatter, appears to hold, over about three orders of magnitude in Li^^ m) from the bright 
fluxes sampled by the IRAS sources (Rush et al. 1993) to the faintest ISO fluxes sampled in 
the HDF-S and HDF-N fields, with more infrared luminous galaxies having on average larger 
Li5^ m /L R ratios. In first approximation, as suggested by the average SEDs shown in Figure 
2, the ratio L 15Atm /L R can be interpreted as an indication of the relative importance between 
bursting and more quiescent emission (see also Rowan- Robinson et al. 2003). In Figure 3 we 
plot the rest-frame Li 5fim /L^ versus Li 5fJim for our galaxies (161 objects), where Li 5fJim and 
L R are the luminosities (uL u ) at 15-/im and in the i?-band, respectively. The majority of 
our galaxies lie between the ratio values found for M51 and M82. On the basis of the plot, 
we have assumed log(Li5 Mm /LR)~— 0.4 as the nominal separation between the normal and 
the starburst galaxy populations (dot-dashed line, see Section 5 for a discussion). This leads 
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to define subsamples of 81 spiral galaxies (80 at z <0.4) and 80 (70 at z <0.4) starburst 
galaxies. We did not use the spectroscopic classification to separate the two populations, 
since dust can significantly depress the emission features, which are mainly used to classify 
starburst galaxies, thus possibly leading to misleading classifications (see Rigopoulou et al. 
2000; Franceschini et al. 2001). This is evident also from Figure 3, where the different 
spectroscopic classes are not clearly segregated in different regions of the plot. In any case, 
it is reassuring that almost all the objects (12/13) spectroscopically classified as e(b) galaxies 
(spectra with very strong emission lines) do indeed have log(L 15Atm /L R ) > — 0.4 and all the 
objects (4/4) classified as k galaxies (elliptical-like spectra) have log(Li 5Mm /L R ) < — 0.4. 

A least square fitting procedure applied to the data leads to the empirical relation: 

log(Li 5 ^ m /L R ) = 0.50 logLi 5 ^ m - 5.4 (1) 

with a dispersion of ~0.28 dex. The best-fit and its la bounds are shown in the figure 
as solid and dotted-lines, respectively. This relation will be discussed in detail in the next 
section. 

In Figure 4 we show the — z diagram. The normal and the starburst galaxies, 

as defined above, are represented by empty and filled circles, respectively. The luminosity 
of each galaxy has been computed using the K-correction appropriate to its classification. 
Considering the galaxies up to z = 0.4, the median redshift of the whole sample (spi- 
ral plus starburst galaxies) is z me d— 0.18, while the median luminosity is Li 5Mm ~lO 9 ' 8 L . 
The majority of the normal galaxies lie at relatively low redshift (z mcd c^0.1A) and luminos- 
ity (L 15/jm ~lO 9 ' 5 L ), while the active population is characterized by higher median values: 
z med ~0.23 and L 15Atm ~lO lo L . 

4. The Mid-IR Luminosity Function 

4.1. Estimator and Selection Effects 

The luminosity function and its evolution have been estimated using the parametric, 
unbinned, maximum likelihood method described in Marshall et al. (1983). We consider 
three different selection effects affecting our data: the 15-/im, the optical i?-band and the 
spectroscopic limits. 

The 15-/im selection effect has been corrected for by weighting each data pair (z,Z/i5 Mm ) 
by its 15-/xm effective area (Q(z, Li^^ m ) = Q(S)). The completeness functions for SI and 
Sl_5 are given in Table 1 of La Franca et al. (2004); for S2, we consider only fluxes greater 
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than 1 mJy with an associated completeness of 93 % for 1 < S < 2 mJy and 100 % for 
S > 2mJy (Pozzi et al. 2003). 

The optical limits, reported in Section 2.2, have been taken into account by introducing 
a function Q(z, Li 5lim ) which represents the probability that a source, with a given 15-/xm 
luminosity Li 5 ^ m and redshift z, had an i?-magnitude within the limits of the sample. To 
calculate the probability of any optical luminosity (Lr) given a specific data pair (z,Li 5lxm ) we 
use the relation between the 15//m and the optical luminosities together with its dispersion, 
assumed to be Gaussian (Eq. 1). The i?-band K-correction described in the previous section 
has been used to relate the optical rest-frame luminosity to the observed magnitude. 

The spectroscopic selection has been considered by weighting each triplet (z, Li 5fJjm , Lr) 
by the spectroscopic completeness in the corresponding optical interval. This correction 
is not significant because of the high spectroscopic completeness in the considered optical 
intervals (>95 % , see Section 2.2). 

Following Marshall et al. (1983) the function to be minimized can be written as S = 
— 21n£, where C is the likelihood function: 

N 

S = -2j2^<f>(zi,Li) (2) 



+ 2 J J cj>(z,L)n(z,L)e(z,L) ^ dz dlogL 



where L is the luminosity at 15 /im (L 15 ^ m ), N is the total number of sources in the three 
samples, Q(z, L) is the available area of the sky for an object with luminosity L at redshift 
z, Q(z, L) is the optical correction factor, (^) the differential volume element and </>(zi, Li) 
is the luminosity function. 

To optimally combine the information from the three different samples (Sl_rest, Sl_5 and 
S2), we follow the formalism described in Avni & Bahcall (1980). Each factor of the double 
integral of Eq. 2 is the sum of three terms: 

n(z,L)Q(z,L) = A sl C sl (z,L)e sl (z,L) (3) 
+ A S i 5 Csi 5 (z,L)Qsi5(z,L) 
+ A S2 C S 2(z,L)Q S2 (z,L) 

where A S i, A S i5 and A S2 are the areas of the three fields (3.55, 0.55 and 0.12 deg 2 for SI, 
Sl_5 and S2 respectively), Csi, Csi5,Cs2 are the three completeness functions and Qsi, ©si5 
and G>s2 are the three optical factors. 

Since we do not have enough data to assess different parametric forms, we decided to 
parameterize the luminosity function <f>(L) using the form suggested by Saunders et al. (1990) 



- 8- 



as a good description of local 60-/zm luminosity function of IRAS galaxies. By using a large 
sample of sources (~2800 objects) Saunders et al. (1990) found that a Schechter function 
was too narrow to describe IR selected sources and a better fit could be achieved using the 
function: 



dN(L,z = 0) 



dz dlogL 

l-a 



(4) 



exp 



2a 2 



In our likelihood analysis, we have searched for the best fitting parameters of the lo- 
cal luminosity function and simultaneously tried to constrain the evolution of the Mid-IR 
galaxies to reproduce the observed distribution of our data in the (z,Li 5fim ) plane. 

As suggested by Franceschini et al. (2001), the shape of the observed source counts 
(Euclidean from IRAS to a few mJy, followed by a sharp upturn at fainter fluxes) favors 
the hypothesis of strong evolution for only a fraction of the whole population, with the 
remaining galaxies giving rise to the Euclidean behavior. We have thus constrained our 
model by assuming that the spiral population does not evolve, while allowing the active 
galaxies (starbursts) to evolve both in density and in luminosity, according to 

<t>(L,z)=g(z)<i>(L/f(z),0) (5) 

parameterizing the two evolutions with two power-laws: g(z) = (l + z) kd and f(z) = (l + z) kl . 

Since the likelihood method determines the shape and evolution of the luminosity func- 
tion, but not the overall normalization, we have normalized the two luminosity functions 
(for normal spiral and active galaxies), by requiring agreement between the predicted and 
the observed total number of sources for each class of objects. 



4.2. Results 

Figures 5a,b show the results of our ML best fit to the luminosity function of spiral 
and starburst populations in two different redshift bins: 0.0 < z < 0.2 (^ mcan ~0.12) and 
0.2 < z < 0.4 (^ mcan ~0.27). The plotted data points correspond to the space densities of the 
observed sources computed independently with the l/V max formalism (Schmidt 1968; Felten 
1976). 
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The best-fitting parameters for the local luminosity functions of the two populations 
are reported in Table 1. The quoted errors correspond to the la confidence limit for each 
parameter calculated while allowing all the other parameters to vary (A£=l, see Lampton 
et al. 1976). While for the starburst population we have allowed the evolution parameters 
to vary, as said before, we have assumed no evolution for the spiral component. For this 
reason, no error is reported for its evolution parameters. Support to this hypothesis is given 
by the V/V max test (Schmidt 1968): in fact, under the hypothesis of no evolution we find 
V/Vmax = 0.55 ± 0.03 for the normal spiral population and V/V max = 0.64 ± 0.03 for the 
starburst population (>4<r evidence of evolution). The latter value becomes 0.52 ± 0.03 
assuming the evolution rates reported in Table 1. 

Because of the relatively small number of objects in each population, the parameters 
derived from our maximum likelihood procedure are not very well constrained. This is true 
in particular for the evolutionary rates of the starburst population. In fact, although our 
data indicate a strong evolution for this population, the uncertainties on the evolutionary 
parameters are large due to the limited redshift interval covered by our surveys. For this 
reason, in order to better test the evolution, we have considered also other observables like 
source counts and redshift distributions at higher z, as will be discussed in Section 6. 

As shown in Figure 5a, the two populations sample different regions of the luminosity- 
redshift plane. In particular, the normal spiral population, which mainly comprises low 
redshift and low luminosity galaxies, samples well the faint end of the luminosity function, 
allowing an accurate determination of the a slope (~20 % uncertainty, see Table 1). On 
the other hand, the starburst population samples well the high luminosity, moderate-high 
redshift regions. This allows the knee of the luminosity function to be better sampled and the 
L* and a parameters to be determined quite accurately. On the other hand, since we have 
not starburst galaxies at low luminosities (see Fig. 5), our data do not allow to constrain, 
for this population, the a slope, which has been fixed to 0.0 as reported in Table 1. 

To test the consistency between our observed z, Li 5fJim distribution and that predicted 
from our parametric model we have performed a 2-dimensional Kolmogorov-Smirnov test 
(2D-KS, see Peacock 1983; Fasano & Franceschini 1987). The 2D-KS test gives >10 % 
probability that the observed data are randomly sampled from the distributions predicted 
by the fitted LF. 
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4.3. The unidentified objects 

In our maximum likelihood procedure, we have considered only objects with R mag 
brighter than our adopted limits and with z<§A (150 objects). We have then corrected 
the remaining redshift incompleteness inside the considered magnitude range by applying 
weights to each galaxy with spectroscopic redshift. In this way, we have assumed that the 
objects with no z (but within the R mag limits) have the same properties (i.e. follow the 
same Li 5fim / Lr vs. Li 5fim relation) as the spectroscopically identified objects with similar R 
magnitude. Since the redshift completeness in the considered magnitude intervals is always 
very high (>95%), we are confident that the uncertainties introduced are negligible. 

A more delicate task is to deal with the 15-/xm sources which have an R magnitude 
fainter than the selection limits or are unidentified (125 objects, ~36 % of the non-stellar 
sample, see La Franca et al. 2004 and Fig. 1). These sources could be either higher redshift 
sources, or optically less luminous galaxies at redshift similar to those of the spectroscopically 
identified sample. In the first case, they are expected to follow the L 15Atm /L^ — L 15/jm relation 
as the other galaxies, while this would not the case for the second hypothesis. 

Following La Franca et al. (2004), the first scenario seems more probable. The con- 
siderations discussed in La Franca et al. (2004) for this choice were mainly two: first the 
LiZfmJ Lr vs. Li 5 ^ m seems to be a valid relation for all the Mid-IR surveys, from the ISOCAM 
ultra-deep to the local IRAS surveys; second, by estimating the redshift of the unidentified 
objects on the basis of the observed log z — R relation, the computed Li 5fim and Lr appear 
well consistent with the assumed relation. In Figure 6, the observed redshift distribution of 
the spectroscopically identified objects and the estimated distribution (on the basis of the 
log z — R relation) of the unidentified sources is compared with the distribution predicted by 
our model (by extrapolating the results to high-z, see Section 6). In Figure 6, the agreement 
between the model predictions and the observed data (including the estimated distribution) 
is quite impressive. In particular, the estimated ^-distribution of the unidentified sources is 
expected to fill exactly the secondary peak predicted by the model. 

A further evidence supporting the hypothesis that most of the unidentified sources 
should belong to the same population as the identified sources, but with higher redshift 
(z = 0.5 — 1.5) is supplied by the photometric redshift technique. The photometric redshifts 
have been estimated and presented in the Final Band-merged ELAIS Catalogue (Rowan- 
Robinson et al. (2003)). The final Band-merged catalogue contains 1636 15-fim sources, 
136 of which with i?>20.5 (optical limit of our larger sample, Sl_rest, see Section 2.2) and 
enough optical data to determine photometric redshifts. The resulting photometric redshift 
distribution for all the ELAIS fields is in agreement with our hypothesis, showing that most 
of the 15- /im sources with associated faint optical galaxies are at significantly higher redshift 
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than sources with brighter optical counterparts (z mean = 0.75 instead of z mean = 0.2 of the 
present optical bright sample). 

The bimodal behavior of the ^-distribution of the 15-/xm sources shown in 6 is caused 
by a combination of different causes: the shape of the LW3 K-correction (see Fig. 2), the 
high evolution rates found for the starburst population and, finally, the slope of the LF for 
high 15-/ZH1 luminosities. The dip in the redshift distribution around z ~0.5 is probably the 
cause of the gap around R ~21 of the distribution of the optical counterparts of the 15-//m 
sources (see Fig. 1 and the discussion in La Franca et al. 2004). 

5. The local luminosity function 

In Figure 7 the local luminosity function (LLF) of 15-fim galaxies (excluding AGNs) 
estimated in this work is compared with other determinations derived at different Mid-IR 
bands and converted to 15-//m using our SEDs (M51 or M82, depending on galaxy type). 

Our estimate of the LLF has been done by extrapolating the result of the maximum 
likelihood method to z = (solid line). To check the over-all normalization we have used the 
l/V max formalism, 'de-evolving' each galaxy according to the derived evolution coefficients 
(filled circles). In Table 2 the 1/V max LLF is listed, with L 15Atm defined as vL v and bin width 
5 log(Li5^ m ) = 0.4. The results of the two methods are well consistent with each other and 
the final V/V max value for the total sample is V A /K naa; =(0.53±0.02). 

The open triangles are an estimate of the LLF of all galaxies (excluding Seyl and Sey2) 
based on the 12-/xm catalogue of Rush et al. (1993). We have computed this LLF using the 
l/Vmax method, selecting all galaxies with Su^m > 300 mJy for which the 12-fim, optical 
and spectroscopic completeness are 100 % (see Rush et al. 1993). The dashed-line is the LLF 
assumed by Franceschini et al. (2001) and is based on the 12-/xm LLF computed by Fang 
et al. (1998) and re-adapted by Xu et al. (1998) in the low luminosity regime. Franceschini 
et al. (2001) assume that starburst galaxies contribute ~10 % of the 12-/im LLF at all 
luminosities. Therefore we have converted their 12-/zm LLF to 15-/im by using the M82 and 
M51 SEDs for ~10 % and ~90% of the 12-/zm LLF, respectively. Moreover, since in the 
Franceschini et al. (2001) model the active population includes also the Sey2 galaxies, we 
have obtained only the starburst contribution by subtracting the LLF of Sey2 computed by 
Rush et al. (1993). Finally, the dot-dashed line is the LLF computed by Xu et al. (2001), 
which is based on the 25-/zm sample of Shupe et al. (1998). This LLF is the sum of an 
actively starforming and a normal population, defined on the basis of the IRAS colours. 

As shown in Figure 7, the three independent estimates are reasonably well consistent 
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with each other. The differences in the low luminosity regime between the Rush et al. 
(1993) and the other determinations is possibly caused by the effect of local inhomogeneities 
(particularly the Virgo super-cluster) in the IRAS survey, as suggested by many authors 
(Fang et al. 1998; Xu et al. 1998). At high luminosity the tendency of the Xu et al. (2001) 
LLF to be higher than the others is probably due to contamination from Sey2 objects, 
showing IR colours similar to those of starbursts. It must be underlined the agreement 
between the overall normalization of our determination (based on the ISO data) and the 
other estimates based on IRAS data. This is indicative of a great accuracy both in the ISO 
calibration achieved with the Lari method (see Lari et al. 2001) and in the completeness 
corrections applied to our data. 

While the total determinations of the Mid-IR LLF for galaxies agree so well, their 
subdivisions into different populations (starburst and normal galaxies) made by different 
authors do not show the same level of consistency. In Figure 8 our LLF for the two galaxy 
populations (thick lines) are compared to those derived by Rush et al. (1993), Franceschini 
et al. (2001) and Xu et al. (2003) (top, middle and bottom panel respectively). The starburst 
populations are shown as dashed lines and the spiral ones as dot-dashed lines. The Rush 
et al. (1993) 12-fim LLF for normal galaxies (shifted to 15-fim through the M51 SED) is 
almost identical to our determination for the same population. Instead, the LLF of their 
liner + starburst component (the latter defined as sources with high FIR luminosity) is 
significantly different from that estimated for our starburst population (basically it is much 
flatter than ours at L 15Atm >10 10 L Q and has a much higher volume density for the highest 
luminosity objects, Li 5 ^ m ~ 10 11 L Q ). The Franceschini et al. (2001) LLF for galaxies are 
very similar to ours, both for starburst and normal spirals (also in this case, the LLF of 
Sey2 computed by Rush et al. (1993) has been subtracted from the active component of 
Franceschini et al. (2001)). The Xu et al. (2003) populations' subdivision, derived from 
IRAS colours of 25-fim selected sources, is totally different, almost opposite to ours. In 
fact, their starburst LLF is higher than the normal galaxy one for L > 10 9 L & and its high 
luminosity slope is rather flat, producing a significant local contribution of starburst galaxies 
even at L > 10 10 5 - 10 11 L & . On the contrary, the Xu et al. (2003) normal galaxy LLF is 
lower and steeper at high luminosities than the starburst one, and is more similar to our 
determination for starbursts. This significantly different subdivision, together with different 
evolutionary schemes proposed, is probably the main cause of the large difference in the 
model predictions between our work and that of Xu et al. (2003) . 

The populations' separation adopted in our work is based on a physical property of 
galaxies, like the ratio between Mid-IR and optical luminosity. The choice of the value of 
the ratio (R = Li 5 ^ m / Lr) chosen to divide the starburst from the normal spiral population 
is mainly based on two reasons. First, as previously said, it is intermediate between the value 
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observed for M51 (normal spiral) and that observed for M82 (starburst galaxy). Second, it is 
the value that best allows to reproduce all the observables, from the normalization of the LLF 
to the redshift distributions at low- and high-z and source counts at all flux levels. A smaller 
fraction of the normal, non-evolving, population (logi?< —0.6), while would have allowed to 
better model the sharp increase of the ELAIS-S1 counts (since a lower contribution of the 
spiral component would have been produced at mJy level), it would have underestimated the 
total local luminosity function. On the contrary, a larger fraction of the normal population 
(logi?> — 0.2) would have enhanced the contribution of the quiescent component in the 
sources counts, predicting a too smooth behavior with respect to the counts shape observed 
in the ELAIS-S1 field. 



6. Discussion on evolution 

As discussed in Section 4.2, the redshift range sampled by the galaxies used to determine 
the LF of the two populations is too narrow (0.0 < z < 0.4) to constrain the evolution 
rates with a high degree of confidence. This is shown in Figure 9, where the contribution 
of the galaxies used in the fitting procedure (dashed line) and of the unidentified objects 
(dot-dashed line) to the total observed source counts (solid line) are illustrated separately. 
The spectroscopic sample contributes only marginally to the large excess with respect to 
the Euclidean expectations at ~ lmJy, which is indeed dominated by optically unidentified 
sources. For this reason, to test our evolution rates, we have compared the predictions of 
our model with other observables, extending the analysis to lower flux densities and higher 
z than those reached by our survey (using the source counts and the HDF-N data-set). 
We have extrapolated the model results to higher redshift, given the evidences reported in 
Section 4.3 that all the 15-/xm sources belong to the same population of galaxies. 



6.1. Comparison between model predictions and data 

We find that the agreement between the model predictions and the observables is very 
good. In particular, the starburst population must evolve with the evolution rates found 
(ki ~ 3.5 and k d ~ 3.8, see Table 1) up to z break ~ 1 and no additional evolution at 
z > zbreak- Of course, the k\ and kd values depend on how the different galaxy populations 
emitting in the Mid-IR band are separated. However, once the subdivision is fixed, the 
evolution rates are well determined, given the large number of observables to be fitted. 

In Figures 10a,b the total (normal spiral+starburst) luminosity functions predicted by 
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our model in the ELAIS southern fields (our sample) and in the HDF-N field, are com- 
pared with the data. For the shallower ELAIS fields, the observed and predicted luminosity 
functions have been computed in two low-redshift bins: 0.0 < z < 0.2 (^ mcan ~0.12) and 
0.2 < z < 0.4 (z mean ~0.27). For the deeper HDF-N, the observed and predicted luminos- 
ity functions are computed in two higher-redshift bins: 0.4 < z < 0.7 (z mean ~0.55) and 
0.7 < z < 1.0 (^ mcan ~0.85). In the latter redshift bins, the model is an extrapolation to 
higher redshift of our best-fitting model; the data points are from Xu et al. (2003) (private 
communication), where they were computed using the V max formalism. 

In Figure 11 the redshift distribution of sources with S > 0.1 mJy predicted by our 
model is superimposed to the data to the same flux density limit observed in the HDF-N. 
The agreement between the observed and the modeled distributions is very good, both in 
shape and normalization, with both distributions showing a peak around z ~ 0.9 — 1.0. 

In Figure 12 the observed and the predicted differential source counts are compared. 
Our model well reproduces the trend observed from the IRAS flux densities down to the 
ultra-deep survey limits. It is in perfect agreement with our ELAIS data at fluxes S i5lJim < 1 
mJy and >3 mJy, while it is slightly higher and smoother than our data in the critical 
interval 1 1 Si 5fMm < 3 mJy, where the counts start diverging from no evolution expectations and 
data from different surveys show the larger differences. Our model is however intermediate 
between, and consistent within the errors with both ELAIS data and those of the other survey 
as reported by Elbaz et al. (1999) and Metcalfe et al. (2003). The sharp upturn shown by 
the ELAIS source counts around 2 mJy could be well reproduced by Gruppioni et al. (2002) 
by introducing a luminosity cut-off in the local luminosity function of starburst galaxies at 
L = 1O 1O ' 8 L . Although our data cannot either confirm or rule out this hypothesis, in the 
present work we have chosen not to introduce any "artificial" constraint in our maximum 
likelihood luminosity function determination. 



6.2. Comparison with other evolutionary models 

The idea of modeling the Mid-IR source counts and luminosity function by dividing the 
sources into different populations, following different evolutionary schemes, was first pro- 
posed by Franceschini et al. (2001). The Franceschini et al. (2001) and our local population 
subdivisions are similar (see Section 5), although the evolutionary rates required by Frances- 
chini et al. (2001) for the starburst population are slightly higher than ours: ~(1 + z) 3 ' 8 in 
luminosity and ~(1 +z) A in density. For this reason, the source counts predicted by Frances- 
chini et al. (2001) are somewhat higher and smoother than ours, especially in the flux density 
interval 1 < S < 5 mJy, and their modeled redshift distribution, although similar at faint 
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flux densities (i.e. Si5 Mm >0.1 mJy in the HDF-N), shows a significant peak around ^~0.8 — 1 
even at relatively high flux densities (>2 mJy and up to 10 mJy), not observed in the data. 
Despite these differences, the Franceschini et al. (2001) model is the model which, among all 
those existing in literature, produces results more similar to ours. 

For example, the recent model of Xu et al. (2003) predicts source counts that are not 
only significantly higher than our expectations, but also higher than all the 15-//m source 
counts derived by the different Mid-IR surveys, over a large flux density interval (1 < 5*15^™ < 10 
mJy). The large discrepancy between the Xu et al. (2003) and our models might be due 
to a combination of several causes, including a different population subdivision of the local 
luminosity function (see Fig. 8) and different AGN contribution. As shown in the previous 
section, the local luminosity function of starburst galaxies, derived by Xu et al. (2003) (see 
Fig. 8) has a very pronounced high luminosity tail in contrast with our determination. 
This tail causes the high predicted counts at few mJys, since local galaxies, characterized by 
-^i5^TO~lO n L and undergoing a moderate-high evolution, at a typical redshift of z^l would 
have Li5 Mm (,2~l)~8 — 10Li 5fim (z = 0) and flux densities in the mJy range. On the contrary, 
local starburst galaxies with typical luminosities of Li 5 ^ m < 1O 9 ~ 1O L (as in our model) would 
have L 15Mrn ~10 11 L Q at z ~ 1 and expected fluxes in the sub-mJy range (S ~0.4 mJy), where 
the bump observed in the differential counts is located. Moreover, the evolution considered 
by Xu et al. (2003) for AGN is higher than that found by I. Matute (private communication). 

Although all the models give rise to different results, it is interesting to note that they 
all agree in the determination of the total local luminosity function (see Section 5), though 
with different population subdivisions and/or evolution hypotheses. 

6.3. Star formation History predicted by model 

The evolving luminosity function model can be used to determine the star formation rate 
density of the Universe. We first calculate the luminosity density at 15-/zm by integrating 
the luminosity function over all the luminosities; then we convert the luminosity density into 
star-formation density by using the star-formation calibrator based on Mid-IR data. To this 
purpose, we have used the calibration given by Mann et al. (2002) for a Salpeter (1955) IMF, 
over the mass range [0.1,100]M Q . Since the Mann et al. (2002) infrared estimator is based on 
the L 60fim bolometric luminosity (SFR(M0yr _1 )=AL A (6O;um)/1.5-lO 36 (W)), L 15fJim has been 
converted to L 60tJim following Mazzei et al. (2001) (LQ 0fim /Li 5lim ~5). The adopted value is 
consistent within ~20-30 % with the values of M51 and M82 (see Fig. 2). In the same way, 
we have also translated into SFR density the l/V max luminosity function results derived in 
two redshift bins (0.0 < z < 0.2 and 0.2 < z < 0.4). 
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In Figure 13 we show a compilation of estimates on the star formation rate density 
as a function of redshift from different indicators taken from Somerville et al. (2001), in a 
Q m = 0.3, Q\ = 0.7 cosmology in units of hM yr~ 1 Mpc~ 3 . The UV data have been corrected 
for dust extinction following Somerville et al. (2001), while the Ha data have been corrected 
for dust extinction by the original authors. The estimates derived from Mid-IR data in the 
HDF-S by Mann et al. (2002) and from radio data by Haarsma et al. (2000) have also been 
added, after conversion to the adopted cosmology. The prediction of our model is shown as 
a solid line and the l/V max data points as filled circles. 

Our model predicts a trend for the star-formation density similar to the results obtained 
in other bands, with a rapid increase from z ~0 to z ~1, followed by a flat plateau at high 
z. The actual measured data points at z < 0.4 are of particular interest, since they provide 
an estimate of the star-formation density at a relatively low redshift, but not so local to 
be affected by clustering and local dishomogeneities. We found star formation density of 
p=(0.025±0.007)hM Q yr- 1 Mpc- 3 and p=(0.043±0.020)hM Q yr~ 1 Mpc~ 3 at the two average 
redshifts z = 0.12 and z = 0.27, respectively. At 0.4<2<1.0, where the ELAIS data are 
highly incomplete, the model has been constrained by other high- z observables in literature. 
We do not extrapolate the predictions of our model at z >1.3, since the LW3 ISO filter 
does not allow to efficiently sample this range of redshift (see LW3 K-correction in the Fig. 
2). Our model prediction is consistent with the estimates derived from UV, optical and 
Mid-IR data up to z < 0.4. In particular, our results are in excellent agreement with the 
data result obtained by Mann et al. (2002) from the Mid-IR survey in the HDF-S (filled 
downward-pointing triangles). At 0A<z<l.0, the model is significantly higher than the 
extinction corrected UV data, suggesting that the extinction corrections applied could be 
underestimated at those redshifts. The estimates derived from radio data by Haarsma et al. 
(2000) (filled squares) are sistematically higher than our model predictions by about a factor 
of two. However, these data might be overestimated, since the authors considered as star- 
forming galaxies (thus contributing to the star-formation density) also all the unidentified 
radio sources, even at flux densities where the fraction of elliptical radio galaxies could still 
be significant (see Gruppioni et al. 1999). 



7. Conclusion 

We have presented the first direct determination of the 15-//m luminosity function and 
its cosmic evolution for galaxies. 

As previously found by other authors, three populations of sources give rise to the 15- 
/im emission: the normal, the starburst and the AGN populations, characterized by different 
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cosmic evolution. In this work we have analysed the galaxy component only (quiescent plus 
actively starforming). The contribution of the AGN component is discussed in a separate 
paper (Matute et al. 2002; Matute et al. in preparation). 

The analysis is based on data from the ELAIS southern fields survey. The sample 
is composed by ~150 galaxies in the redshift interval 0.0<z<0.4 and covers a large flux 
density range intermediate between the IRAS and the deep ISOCAM surveys (0.5< S < 50 
mJy). Differently from other authors, we have adopted in this work the L 15Mm /LR ratio as 
a criterion to separate the quiescent, non-evolving and the starburst, evolving populations. 
This criterion, suggested by the existing correlation between Li5^ m /LR and the amount of 
activity in galaxies, is a posteriori supported by the results of the V/V ma x analysis on the 
two populations defined on this basis. 

The main results of our analysis are: 

1. In the ML analysis we have simultaneously fitted both the evolution rates and the 
shape parameter of the local LF for both the spiral and the starburst populations. We 
have assumed that the spiral population does not evolve, while we have let to evolve 
the starburst population both in luminosity and in density. Since the two populations 
sample different luminosity ranges, we have obtained an accurate determination for the 
faint end of the LLF for the quiescent component, while the knee and the a parameters 
of the LLF are better constrained for the starburst one. The evolution found for the 
active population is ~ (1 + z) 3,5 in luminosity and ~ (1 + z) 3 - 8 in density, up to 

Zbreak ~ 1- 

2. Our total 15-fim LLF is in agreement with previous determinations derived from the 
IRAS data. On the contrary, the LLFs for different populations derived by different 
authors have not the same level of consistency. While in our subdivision the quiescent 
population is expected to dominate locally over all the luminosity range, in other mod- 
els (i.e. Xu et al. 2003) the starburst population dominates locally at high luminosities, 
leading to a large discrepancy in the model predictions. 

3. To test the evolution parameters with higher degree of confidence, we have compared 
our model predictions with all the observables existing in literature, over all the z 
and flux ranges (source counts, luminosity functions, ^-distributions). Our best-fitting 
model well reproduces all the observables. In the critical interval 1<S<3 mJy, where 
the source counts from different surveys show the larger discrepancies, our model is 
intermediate between the data from ELAIS-S1 (Gruppioni et al. 2002) and the data 
from the deep surveys (Elbaz et al. 1999). On the other hand, in the flux range 
3tS< 10 mJy, our model is well consistent with existing data, differently from the Xu 
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et al. (2003) model, whose predicted differential sources counts are at least a factor of 
3 higher than the data. 

4. Using the evolutionary model found for the 15-//m galaxies and the data points from the 
l/V max LF analysis, we have estimated the star-formation rate density. The redshift 
range sampled by our data (0.0 < z < 0.4) is of particular interest, since it provides an 
estimate of the star-formation at relatively low redshift, but not so local to be affected 
by clustering and local dishomogeneities. We find p=(0.025±0.007) hM Q yr _1 Mpc -3 
and p=(0.043±0.020) hM yr _1 Mpc _3 at the two mean redshifts z = 0.12 and z = 0.27, 
respectively. At ztOA our model predictions are well consistent with other estimates 
derived from UV, optical and Mid-IR data. At higher redshift our model predictions 
are significantly higher than the UV extinction corrected data and lower by about a 
factor of two than the estimates derived from radio data by Haarsma et al. (2000). 



F.P. was partially supported by the European Program HUMAN POTENTIAL (Con- 
tract Number HPMT-CP-2000-00096). F.P. would like to thank Lucia Pozzetti for helpful 
discussion. 

REFERENCES 

Aussel, FL, Elbaz, D., Cesarsky, C. J., & Starck, J. L. 1999, in ESA SP-427: The Universe 
as Seen by ISO, 1023 

Avni, Y., & Bahcall, J. N. 1980, ApJ, 235, 694 

Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 

Dressier, A., Smail, I., Poggianti, B. M., Butcher, H., Couch, W. J., Ellis, R. S., & Oemler, 
A. J. 1999, ApJS, 122, 51 

Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H., Franceschini, A., Fadda, D., & Chary, 
R. R. 2002, A&A, 384, 848 

Elbaz, D., et al. 1999, A&A, 351, L37 

Forster Schreiber, N. M., Sauvage, M., Charmandaris, V., Laurent, O., Gallais, P., Mirabel, 
I. F., & Vigroux, L. 2003, A&A, 399, 833 



-19- 

Fang, F., Shupe, D. L., Xu, C, & Hacking, P. B. 1998, ApJ, 500, 693 
Fasano, G., & Franceschini, A. 1987, MNRAS, 225, 155 
Felten, J. E. 1976, ApJ, 207, 700 
Flores, H., et al. 1999, ApJ, 517, 148 

Franceschini, A., Aussel, H., Cesarsky, C. J., Elbaz, D., & Fadda, D. 2001, A&A, 378, 1 
Franceschini, A., et al. 2003, A&A, 403, 501 

Gronwall, C. 1998, in ASP Conf. Ser. 146: The Young Universe: Galaxy Formation and 
Evolution at Intermediate and High Redshift, 96 

Gruppioni, C, et al. 1999, MNRAS, 305, 297 

Gruppioni, C, Lari, C, Pozzi, F., Zamorani, G., Franceschini, A., Oliver, S., Rowan- 
Robinson, M., & Serjeant, S. 2002, MNRAS, 335, 831 

Gruppioni, C, Mignoli, M., & Zamorani, G. 1999, MNRAS, 304, 199 

Haarsma, D. B., Partridge, R. B., Windhorst, R. A., & Richards, E. A. 2000, ApJ, 544, 641 

Hauser, M. G., et al. 1998, ApJ, 508, 25 

Hughes, D. H., et al. 1998, Nature, 394, 241 

La Franca, F., et al. 2004, AJ, in press (astro-ph/0403211) 

Lagache, G., Abergel, A., Boulanger, F., Desert, F. X., & Puget, J.-L. 1999, A&A, 344, 322 
Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177 
Lari, C, et al. 2001, MNRAS, 325, 1173 
Mann, R. G., et al. 2002, MNRAS, 332, 549 

Marshall, H. L., Tananbaum, H., Avni, Y., & Zamorani, G. 1983, ApJ, 269, 35 
Matute, I., et al. 2002, MNRAS, 332, Lll 

Mazzei, P., Aussel, H., Xu, C, Salvo, M., De Zotti, G., & Franceschini, A. 2001, New 
Astronomy, 6, 265 

Metcalfe, L., et al. 2003, A&A, 407, 791 



-20- 



Oliver, S., et al. 2000, MNRAS, 316, 749 
Peacock, J. A. 1983, MNRAS, 202, 615 
Pozzi, F., et al. 2003, MNRAS, 343, 1348 

Puget, J.-L., Abergel, A., Bernard, J. -P., Boulanger, F., Burton, W. B., Desert, F.-X., & 
Hartmann, D. 1996, A&A, 308, L5 

Rigopoulou, D., et al. 2000, ApJ, 537, L85 

Rowan-Robinson, M. 2001, ApJ, 549, 745 

Rowan-Robinson, M., & Crawford, J. 1989, MNRAS, 238, 523 

Rowan-Robinson, M., et al. 2003, preprint (astro-ph/0308283) 

Rush, B., Malkan, M. A., & Spinoglio, L. 1993, ApJS, 89, 1 

Salpeter, E. E. 1955, Vistas in Astronomy, 1, 283 

Saunders, W., Rowan- Robinson, M., Lawrence, A., Efstathiou, G., Kaiser, N., Ellis, R. S., 
& Frenk, C. S. 1990, MNRAS, 242, 318 

Schmidt, M. 1968, ApJ, 151, 393 

Shupe, D. L., Fang, F., Hacking, P. B., & Huchra, J. P. 1998, ApJ, 501, 597 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 103 

Somerville, R. S., Primack, J. R., & Faber, S. M. 2001, MNRAS, 320, 504 

Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., & Pettini, M. 1999, ApJ, 
519, 1 

Tresse, L., & Maddox, S. J. 1998, ApJ, 495, 691 

Treyer, M. A., Ellis, R. S., Milliard, B., Donas, J., & Bridges, T. J. 1998, MNRAS, 300, 303 
Xu, C, et al. 1998, ApJ, 508, 576 

Xu, C, Lonsdale, C. J., Shupe, D. L., O'Linger, J., & Masci, F. 2001, ApJ, 562, 179 

Xu, C. K., Lonsdale, C. J., Shupe, D. L., Franceschini, A., Martin, C, & Schiminovich, D. 
2003, ApJ, 587, 90 



This preprint was prepared with the AAS IATgX macros v5.2. 



-21 - 



Table 1: 15-//m Luminosity Function parameters from the ML analysis. 



Population 


a 


a 




log 0* 


k 


fed 


normal spirals 


i - lu -0.25 


U.O_ 2 


o O+0.7 
O-O-0.9 


-2.45 


0.0 


0.0 


starbursts 


0.0 (fixed) 


n qq+0.025 

u.oy_ .025 


o o+0.3 
O-O-0.2 


-3.53 




q o+2.0 
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Table 2: 15-//m Local luminosity Function from the l/V ma x analysis. 



\og[vL u /L & ] 


log [(j)(Mpc 3 mag 1 


la Error 


7.8 


-1.94 


(+0.52)(-0.76) 


8.2 


-2.27 


(+0.24) (-0.64) 


8.6 


-2.37 


(+0.16)(-0.25) 


9.0 


-2.84 


(+0.12)(-0.17) 


9.4 


-2.81 


(+0.08)(-0.09) 


9.8 


-3.24 


(+0.06)(-0.08) 


10.2 


-4.25 


(+0.11)(-0.14) 


10.6 


-5.21 


(+0.23)(-0.53) 
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Fig. 1. — i?-band magnitude as a function of 15/xm flux density for all 15/im sources in 
the ELAIS southern fields. Sources with redshift are shown as filled symbols, sources with 
likely optical counterparts but without redshift are shown as empty symbols and sources 
without optical counterparts are shown as lower limits (arrows). Circles stand for sources 
from Sl_rest; triangles for sources from Sl_5 and squares for sources from S2. The three 
lines represent the three magnitude thresholds considered in this work (dotted-line: Sl_rest; 
dashed line: Sl_5; dot-dashed line: S2). 
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Fig. 2. — Top: Spectral energy distributions (SED) adopted for our galaxies. The dashed-line 
corresponds to M51, a low luminosity inactive spiral galaxy, while the solid-line corresponds 
to M82, considered as the star-forming galaxy prototype. The SEDs have been taken from 
the GRASIL code (Silva et al. (1998)). In the range from 5 to 18 ^m the SED for M82 
is the observed ISOCAM CVF spectrum. The vertical dot-dashed lines correspond to the 
i?-band and LW3 filter trasmissions at the half-maximum-transmission values. Bottom: The 
K-correction as a function redshift in the LW3 and i?-band filters for M51 (dashed-lines) and 
M82 (solid-lines). 
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Fig. 3. — The Li 5fJim /Ln ratio as a function of Li 5fJim for galaxies in the spectroscopic sample. 
Different symbols correspond to the different spectroscopic classes (see La Franca et al. 2004). 
The large diagonal crosses represent the values for M51 and M82. 



O 



8 



0.0 0.2 



0.4 
z 



0.6 0.8 



Fig. 4. — Rest-frame 15-//m luminosity as a function of redshift for galaxies in the spec- 
troscopic sample. Empty circles represent objects with log(Li 5/wn /LR.) < —0.4 (defined as 
normal galaxies in this work); filled circles represent objects with \og(Li 5fim / Lr) > —0.4 
(defined as starburst galaxies in this work). 
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Fig. 5. — Top: The rest-frame 15-/xm luminosity functions for the two galaxy classes from 
our survey in the two redshift bins 0.0 < z < 0.2 and 0.2 < z < 0.4. The empty circles 
represent the l/V max determination in the redshift interval 0.0 - 0.2, while the filled circles 
are in the interval 0.2 - 0.4. upper plot normal galaxy population; lower plot: starburst 
population. Bottom: Comparison between the observed (dot-dashed histogram) and the 
predicted redshift distributions (solid-line) for the two galaxy classes, upper plot: normal 
galaxy population; lower plot: starburst population. 
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Fig. 6. — Comparison between the observed redshift distribution for the spectroscopically 
identified objects in the total Sl_rest+Sl_5+S2 sample (long-dashed histogram) and the 
model predictions (solid-line). The contributions of the normal spiral and starburst com- 
ponents are shown as dot-dashed and dashed lines respectively. The estimated distribution 
of the unidentified objects as given in La Franca et al. (2004) is also shown as dot-dot-dot- 
dashed line. 
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Fig. 7. — Local 15- jum luminosity function of galaxies. The filled circles are the l/V max 
estimates from this work, while the solid line is the extrapolation to z = of the ML results. 
The dashed line is the 12-/xm LLF from Franceschini et al. (2001) (the Sey2 contribution has 
been subtracted using the LLF for Sey2 computed by Rush et al. 1993). The dot-dashed 
line is the 15-/xm conversion of the LLF computed by Xu at 25-/im. Empty triangles are 
an estimate at 12-/im of the LLF of galaxies, based on the Rush et al. (1993) catalogue. 
The LLF computed at Mid-IR bands different from 15-/xm have been converted to 15-/im as 
explained in the text. 
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Fig. 8. — Comparison between the LLF subdivision into starburst and normal galaxy popu- 
lations from this work and other analyses. The starburst galaxies are shown as dashed lines 
while the spirals as dot-dashed lines. The determinations from this work are thicker. Top: 
comparison with Rush et al. (1993); Middle: comparison with Franceschini et al. (2001) (the 
LLF computed by Rush et al. (1993) for the Sey2 galaxies has been subtracted from the 
active component). Bottom: comparison with Xu et al. (2003). 
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Fig. 9. — The contribution to the total integral source counts (solid line) from galaxies in 
the spectroscopic sample (dashed-line) and from unidentified sources (dot-dashed line) in 
the ELAIS-S1 survey. The areas filled with horizontal and diagonal lines represent the 
68% confidence regions. The confidence region for the total counts have not been shown for 
clarity. 
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Fig. 10. — Top: The rest-frame 15-//m luminosity functions for the total sample of galaxies 
in the southern ELAIS fields (S1+S2). Dashed, dashed-dotted and solid lines represent 
the starburst, normal galaxies and total LFs from our model, upper plot: 0.0 < z < 0.2 
(<2mcan— 0.12). lower plot 0.2 < z < 0.4 (z mean ~0.27) . Bottom: The predicted rest-frame 15- 
fim luminosity functions extrapolated to higher redshift intervals and compared with data 
from the HDF-N survey. The points have been derived by Xu et al. (2003) from 1/V max 
analysis (private communication), upper plot: 0.4 < z < 0.7 (^ mean ~0.55). lower plot: 
0.7 < z < 1.0 (z mean ~0.85). 
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Fig. 11. — Comparison between the observed (dot-dashed histogram) and the predicted 
redshift distributions (solid-line) in the HDF-N field. Sources with S f >0.10 mJy have been 
considered (41 objects; H. Aussel and S. Berta, 2003, private communication). 
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Fig. 12. — Comparison of model prediction and observed 15-/xm source counts. Data points: 
the ELAIS-S1 counts reported by Gruppioni et al. (2002) are plotted as filled stars. The 
new A2390 source counts from Metcalfe et al. (2003) are plotted as open diamonds. Other 
data points from Elbaz et al. (1999): ISO HDF-N (open circles), ISO HDF-S (filled cir- 
cles), Marano Firback (open squares), Marano Ultra-Deep (diagonal crosses), Marano Deep 
(asterisk), Lockman Deep (open triangles), Lockman Shallow (filled triangles). Dashed and 
dashed-dotted lines represent the contributions from starburst and normal galaxies computed 
with our model. 
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Fig. 13. — A compilation of constraints on the star formation history of the Universe mainly 
taken from Somerville et al. (2001). The UV data have been corrected for dust extinction 
following Somerville et al. (2001), while the Ha data have been corrected by the original 
authors. In all cases a Q m = 0.3, Q/y = 0.7 cosmology and a Salpeter (1955) IMF, over the 
mass range [O.1,1OO]M , were assumed. Units are hM yr _1 Mpc -3 . The filled circles and 
the solid curve are from this work (1/V max analysis of ELAIS data and prediction of our 
maximum likelihood model, the dashed and the dot-dashed lines being the contribution of 
starbursts and spirals respectively). The symbols are as follows: empty upward-pointing 
triangles - Gronwall (1998); empty circles - Tresse & Maddox (1998); empty diamonds - 
Treyer et al. (1998); filled downward-pointing triangles - Mann et al. (2002); filled upward- 
pointing triangles - Flores et al. (1999); filled squares - Haarsma et al. (2000); empty squares 

- Cowie et al. (1996); empty downward-pointing triangles - Steidel et al. (1999); filled stars 

- Hughes et al. (1998). 



